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Ion-molecule  reactions,  ion  clusters,  nucleation  phenomena, 
structure  and  bonding  of  small  complexes,  aerosols,  formation 
of  surfaces,  small  particles 


Durinq  recent  years/there  has  been  an  enormous  growth  of  interest 
in  the  properties  of  small  particles  and  in  the  mechanisms  by 
which  they  are  formed.  This  has  been  due  partly  to  the  realiza¬ 
tion  that  findings  in  this  field  can  contribute  to  an  understand¬ 
ing  of  a  very  wide  range  of  fundamental  problems  concerned  with 
phase  transitions,  the  formation,  structure  and  reactivity  of 
surfaces,  the  nature  and  stability  of  the  aerosol  stpte,  and  even 
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20.  Abstract  (continued) 

the  kinetics  of  association  and  unimolecular  reactions,  and  sol¬ 
vation  phenomenon.^  Further  impetus  has  also  arisen  because  of 
the  larqer  number  of  applied  areas  to  which  the  results  have  a 
direct  bearing  including  the  atmospheric  sciences,  obscuration, 
combustion  research,  and  catalysis  to  mention  just  a  few.  In  all 
but  one  of  these  areas  the  system  of  interest  is  that  comprised  of 
particles  dispersed  in  the  gas  phase  commonly  termed  an  aerosol. 

~,sTho  research  proqram  undertaken  in  our  laboratory  was  directed  to 
a  study  of  the  fundamentals  of  nucleation  about  negative  ions, 
with  particular  attention  to  an  elucidation  of  the  structure  of 
small  clusters,  and  the  role  of  chemical  bondinq  in  establishing 
their  stability.  t'sing  a  new  high  pressure  mass  spect romet r ic 
technique,  a  systematic  investigation  of  the  bond  enerqies  and 
entropies  of  numerous  ligand-negative  ion  complexes  were  deter¬ 
mined.  Concurrent  theoretical  work  on  the  bondinq  of  small  com¬ 
plexes  and  ion-induced  nucleation  was  also  undertaken.  The  re¬ 
sults  have  led  to  a  more  complete  understanding  of  nucleation 
phenomena,  and  especially  of  the  importance  of  cluster  structure 
in  the  development  of  surfaces  and  in  effecting  phase  transforma¬ 
tions.  In  addition,  the  work  has  led  to  a  clarification  of  some 
recent  atmospheric  observations  reported  in  the  literature. 
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Summary  of  Findings 

Support  from  the  U.  S.  Army  Research  Office  commenced 
on  10  June  1976  for  a  grant  entitled  "Charged  Atmospheric  Conglom¬ 
erates:  Role  of  Negative  Ions  in  Heteromolecular  Nucleation;  the 

Properties  of  Small  Clusters  and  Atmospheric  Ions."  The  early 
period  of  the  grant  was  devoted  to  reconstructing,  with  numerous 
modifications,  an  apparatus  moved  to  the  University  of  Colorado 
when  the  principal  investigator  accepted  a  position  there.  Fol- 
owing  necessary  testing  procedures,  an  experimental  program 
devoted  to  providing  the  molecular  basis  for  understanding  the 
stability  of  negative  ion  clusters  was  begun.  Related  theoretical 
work  on  ion-cluster  structure,  and  for  assessing  the  validity  of 
conventional  treatments  of  heteromolecular  nucleation, was  also 
done  concomitant ly. 

The  results  of  this  research  have  led  to  the  publication 
of  eight  papers  (currently  in  print  or  in  press)  attached  as 
Appendices  A-H,  presentations  at  39  meetings  and  seminars,  and 
results  for  six  other  journal  articles  which  are  now  in  prepara¬ 
tion.  Since  the  results  are  detailed  in  the  attached  papers,  only 
a  very  brief  summary  of  the  findings  is  given  here.  In  accordance 
with  the  original  proposal,  the  work  was  directed  largely  to 
systems  of  interest  in  the  upper  atmosphere. 

1 .  Studies  of  the  Bonding  of  Negative  Ion  Clusters 

It  is  well  known  that  both  natural  events,  such  as 


solar  flares,  as  well  as  manmade  ones  such  as  a  nuclear  burst,  can 
cause  atmospheric  electron  densities  and  levels  of  ionization  to  be 
increased  to  values  several  orders  of  magnitude  higher  than  those 
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of  the  unperturbed  atmosphere.  Since  the  propagation  of  directed 
electromagnetic  radiation  (i.e.,  radar,  radio  signal,  and  optical 
signals)  is  influenced  by  the  intensity  of  these  levels,  a  deter¬ 
mination  of  the  chemical  and  physical  processes  controlling  the 
rate  of  return  to  normal  atmospheric  conditions  was  deemed  of 
considerable  importance  so  that  the  effect  of  these  perturbations 
on  communication  could  be  properly  assessed. 

Although  considerable  attention  had  been  given  to  a 
study  of  the  homogeneous  chemical  reactions  and  physical  charge 
exchange  processes  which  might  operate  in  the  atmosphere,  by  com¬ 
parison  very  little  was  known  about  the  processes  associated  with 
clusters,  conglomerates,  and  small  aerosols.  Of  particular  inter¬ 
est  is  an  understanding  of  the  properties  and  fate  of  negative  ions. 
Mass  spectrometer  measurements  made  in  the  upper  atmosphere  by 
Krankowsky,  et  al.*  *  '  and  Narcisi^  ^  have  showed  the  existence 
of  an  abundance  of  various  negative  ion  clusters.  The  findings 
of  Krankowsky  and  coworkers  showed  CO^  to  be  a  major  ion,  followed 
by  HCO^  and  NO^j  the  data  also  indicated  the  existence  of  hydra¬ 
tion  complexes  of  these,  as  well  as  the  presence  of  CO^,  and 
02  (H20) ^  *  and  perhaps  Cl  and  NO'^HNO^*  Somewhat  in  conflict  were 
the  results  of  Narcisi  who  reported  the  tentative  identification 
of  N03(H20)n  (n=0-8),  C03(H20)2,  and  C0~.  The  latter  believes  that 
Cl  and  N02  might  originate  from  rocket  contamination,  but  this 
suggestion  is  disputed  by  Krankowsky  and  Arnold. 

Tn  order  to  provide  the  data  necessary  for  interpreting 
the  stability  of  these  complexes,  an  investigation  of  the  thermody¬ 
namic  properties  of  COj  and  HCO^  complexes  with  water  was  first  under¬ 
taken.  Specifically,  extensive  studies  were  made  of  the  reaction: 
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C03(H20)n-l+H2°  «="=*  C03(H20)n 
with  n  ranginq  from  1  to  3,  and 

HCO~  (H20)n_1-*-H20  >  HC0j(H:0)n 

with  n  ranging  from  1  to  4  .  Experiments  were  performed  over  a 
wide  range  of  pressures,  temperatures,  and  electrical  field  ener¬ 
gies  and  it  was  established  that  equilibrium  was  indeed  attained 
in  the  reaction  cell.  The  technique,  experimental  details,  and 
results  are  detailed  in  Appendix  E.  From  measurements  of  the 
temperature  dependencies  of  the  equilibrium  constants,  both  en¬ 
thalpies  and  entropies  were  derived  for  the  reactions.  The 

results  for  reaction  (1)  are  -AH°  ,  *14.1,  13.6,  and  13.1  kcal/mole 

n-i  9  n 

and  -AS^  .  -25.2,  29.6,  and  32.5  e.u.  for  n-1  to  3,  respectively. 

n"l » n 

For  reaction  (2)  the  respective  enthalpy  values  are  -15.7,  -14.9, 
-13.6,  and  -13.4  kcal/mole;  the  corresponding  entropies  are  -24.1, 
-29.1,  -30.2,  and  -33.3  e.u. 

An  interesting  finding  is  that  the  thermodynamic  proper¬ 
ties  of  the  C0~  hydrates  nearly  merge  with  those  of  HCO^  of  one 
less  cluster  size,  beyond  that  of  the  first  cluster  of  CO^.  Of 
particular  importance  is  the  fact  that  the  enthalpies  of  hydration 
do  not  readily  approach,  asymptotically,  the  heat  of  condensation 
for  H20  as  is  often  the  case  for  positive  ions.  (See  Appendix  E, 
Figure  5).  Therefore,  the  hydration  of  these  negative  ions  is 
stronger  than  those  of  many  analogous  simple  positive  ions  which 
do  display  this  tendency.  This  gives  credence  to  the  rocket 
borne  mass  spectrometric  findings  of  Arnold  and  Krankowsky  regard¬ 
ing  the  importance  of  CO~  and  HCOj  hydrates  in  the  mesosphere. 
Details  are  discussed  in  Appendix  F.  As  an  aid  in  understanding 
the  bonding  resulting  from  hydration  of  these  ions,  approximate 


quantum-mechanical  molecular  orbital  calculations  were  performod 
with  the  CNDO/2  method  for  one,  two,  and  three  water  molecules 
bound  to  HCOj.  Unlike  the  situation  for  positive  ions,  hydration 
complexes  of  neqative  ions  develop  through  bonding  of  a  hydrogen 
atom  to  the  ion. 

Intense  interest  also  exists  in  the  ions  N02  and  N03« 
and  we  have  completed  studies  of  their  interaction  with  H20,  as 
well  as  N02  with  HNOj,  CC>2  and  S02#  and  N0~  with  HNC>3.  Results 
of  N0x  bonded  with  H20  and  HNOx  are  discussed  in  Appendix  G.  (See 
Figures  1,  2,  3  and  4  (in  the  Appendix)  for  the  equilibrium 
constants  determined  at  a  wide  variety  of  temperatures. 

Recent  data  acquired  in  our  laboratory  for  the  bond 
energies  and  entropies  of  the  various  negative  ion  cluster  reactions 
are  given  in  Table  I,  along  with  data  for  the  hydration  complexes 
of  CO^  and  HCO^  for  comparison.  The  importance  of  the  nature  of 
the  ligand  on  its  ability  to  bond  to  negative  ions,  and  partic¬ 
ularly  with  regard  to  establishing  a  stable  cluster  species, 
is  seen  by  comparing  data  for  a  number  of  different  ligands 
bonded  to  one  ion.  Referring  to  Figure  1,  equilibrium  constants 
for  a  number  of  ligands  bonded  to  N02  are  seen  to  be  greatly 
different . 

Significant  differences  are  also  observed  in  the  case 
of  various  ions  bound  to  a  single  ligand.  An  example  is  given 
in  Figure  2  which  is  a  plot,  versus  reciprocal  temperature,  of 
the  logarithm  of  the  equilibrum  constants  ci  Cl",  N0~,  S02  and 
S03  bonded  with  SC>2.  On  this  basis,  complexes  of  S02  with  C£  and 
N02  should  be  especially  important  in  regions  where  these  parent 
ions  exist;  however,  complexes  with  S0~  are  unlikely  to  be 


TABLE  I 


THERMODYNAMIC  PROPERTIES  OF  NEGATIVE  ION  CLUSTERS 
RECENTLY  MEASURED  IN  OUR  LABORATORY 


~~  I 

A  Bn  +  B  ^  A  Bn+i 


|  “AH°nf  n>,  (kcal/mole) 

-AS’ 

*n,  n+i 

(gibb/mole) 

Ligand 

ion 

(0,1) 

(1,2) 

(2,3) 

(3,4) 

(0,1) 

(1,2) 

(2,3) 

(3,4) 

c  r 

14.9 

12.6 

11.5 

10.9 

19.7 

20.5 

— 

22.4 

24.8 

r 

II. 1 

9.9 

9.3 

19.3 

20.3 

21.0 

C0'3 

14.1 

13.6 

13.1 

25.2 

29.6 

32.5 

H20 

hco_3 

15.7 

14.9 

13.6 

13.4 

24.1 

29.1 

30.2 

33.3 

no2 

15.2 

13.6 

11.7 

11.6 

23.8 

26.4 

25.8 

29.0 

NO3 

_ 

14.6 

14.3 

■ 

13.8 

25.0 

30.3 

33.2 

cr 

8.0 

19.6 

1* 

5.6 

18.2 

C02 

co-3 

7.1 

21.8 

N0£ 

9.3 

24.2 

SO3 

6.5 

20.7 

cr 

21.8 

12.3 

10.0 

8.6 

23.2 

22.7 

23.1 

23.2 

1 

no2 

259 

9.0 

6.6 

36.8 

16.8 

13.4 

so2 

so2 

24.0 

33.8 

1 

S03 

13.3 

1 

18.9 

I' 

12.9 

10. 1 

9.2 

20.2 

21.6 

24.7 
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important.  Equilibrium  constant*  for  othar  negative  ion  complexes 
with  the  non-polar  molecule  CO^  are  ahown  in  Figure  3.  Apparently, 

the  dipole  moment  of  the  clustering  molecule  plays  a  far  less 

* 

significant  role  in  stabilizing  clusters  of  negative  ions  than  for 
positive  ones. 

The  molecules  HN02  and  HNO^  are  important  in  the  atmos¬ 
pheric  NO^  cycle,  and  interest  exists  in  their  bonding  to  N0~  and 
NO j .  Our  findings  concerning  the  stabilities  of  the  complexes 
also  provide  an  explanation  for  certain  other  unidentified  masses 
detected  by  rocket  borne  mass  spectrometric  experiments. 


Investigations  Pertaining  to  Nucleation  and  Aerosol  Formation 


Heteroqeneous  nucleation,  including  the  subclass  referred 
to  as  heteromolecular  nucleation,  is  the  one  responsible  for  the 
formation  of  atmospheric  aerosols  from  both  natural  and  pollutant 
molecules,  as  well  as  clouds  in  the  lower  and  upper  atmosphere. 

Even  the  formation  of  noctilucent  clouds  in  the  region  of  the  meso¬ 
sphere  is  believed  to  result  from  the  heteromolecular  phenomenon; 
nucleation  about  ions  is  almost  certainly  the  process  responsible 
for  their  formation  ^  .  For  this  reason,  studies  of  ion  induced 
nucleation  provide  direct  information  on  processes  of  atmospheric 
importance,  and  this  has  been  a  major  object  of  the  work  carried 
out  to  date.  By  carefully  selecting  the  ions  and  ligands,  and 
conducting  the  experiments  with  systems  which  can  be  most  readily 
treated  theoretically,  e  great  deal  of  information  on  the  basic 
phenomenon  has  been  obtained.  This  has  led  to  an  understanding  of 
the  limitations  of  the  conventional  theoretical  methods  often 
used  to  treat  nucleation  and  aerosol  formation.  Two  reviews. 
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recently  prepared  by  the  principal  investigator  and  attached  as 
Appendices  A  and  C  detail  the  importance  of  these  studies  in  under¬ 
standing  the  nature  of  the  prenucleation  embryos  to  heteromolecular 
nucleation.  The  data  have  also  found  application  in  interpreting 
the  nature  of  the  interactions  between  ions  and  ligands  such  as 
those  which  occur  in  a  microdroplet,  or  the  solution  phase  consti¬ 
tuting  the  liquid  part  of  an  aerosol  particle,  and  are  thereby 
related  to  solvation  phenomenon. 

As  part  of  the  ARO  sponsored  program  theoretical  con¬ 
sideration  was  given  to  the  applicability  of  classical  theories  of 
homogeneous  nucleation  to  treating  chemically  interacting  systems. 
Experimental  data  obtained  as  part  of  the  ARO  proqram  dealing  with 
negative  ions,  as  well  as  results  from  another  program  underway  in 
our  laboratory,  have  been  useful  in  this  regard.  In  order  to 
formulate  the  rate  of  heteromolecular  nucleation,  it  has  become 
common  to  employ  the  Thomson  equation  which  represents  a  minor 
extension  of  the  Kelvin  relationship  commonly  employed  in  homo¬ 
geneous  theory.  Macroscopic  experiments  in  cloud  chambers  some¬ 
times  confirm  the  utility  of  these  relationships,  but  oftentimes 
suggest  their  inapplicability  in  treating  other  systems.  No 
general  agreement  as  to  validity  of  the  theories  has  been  obtained 
from  investigations  due  to  lack  of  knowledge  concerning  reasons 
for  its  apparent  usefulness  in  some  cases  and  not  in  others. 

The  Thomson  equation  contains  three  terms:  the  first 
accounts  for  the  free  energy  of  vaporization  of  the  condensing 
molecules,  the  second  represents  the  work  done  in  changing  the 


surface  area  upon  condensation  of  additional  molecules  about  an 
ion  cluster,  and  the  third  accounts  for  the  change  in  field  energy 
due  to  the  condensation  of  a  dielectric  film  (constant  c)  about 
the  ion.  All  of  the  capillarity  approximations  implicit  in  the 
use  of  the  Kelvin  equation  for  homogeneous  nucleation  are  involved 
in  the  use  of  the  Thomson  equation;  the  nature  of  an  ion  is  ex¬ 
pressed  only  by  its  ionic  radii.  Consideration  of  the  range  of 
validity  of  the  equation  is  instructive  in  developing  an  under¬ 
standing  of  nucleation  about  ions  and  in  deciding  what  terms  must 
be  sought  in  quantitatively  accounting  for  the  differences  which 
various  ions  display  in  their  ability  to  promote  nucleation. 

An  expression  for  the  free  energy  change  per  molecule, 
^Gn-1  n'  was  reac*ilY  developed  by  taking  the  derivative  of  the 
Thomson  equation  with  respect  to  radius;  derivative  with  respect 
to  temperature  enables  an  evaluation  of  entropy.  Enthalpy  (propor¬ 
tional  to  bond  enercy)  is  then  obtained  by  substitution  of  these 
values  into  well  known  thermodynamic  relationships.  (Refer  to 
Equations  3  and  11,  Appendix  B) . 

It  is  recognized  that  nucleation  will  not  occur  at 
supersaturations  below  100%,  and  this  is  borne  out  by  our  calcula¬ 
tions  made  for  an  ion-water  system  and  shown  in  Figure  1.  (See 
Appendix  B)  Here,  the  energy  barrier  to  nucleation,  A$,  is 
plotted  versus  cluster  number  n  for  varyinq  ratios  of  in  p/p— , 
where  p  is  the  partial  pressure  of  the  condensing  vapor  and  p^ 
the  pressure  above  a  surface  with  an  infinite  radius  of  curvature. 
The  closed  triangular  points  represent  a  saturated  system  where  the 
energy  barrier  is  effectively  infinite;  the  minimum  at  cluster 
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size  10  represents  the  most  abundant  ion-cluster  at  this  concentra¬ 
tion  of  condensable  vapor.  At  higher  pressure  ratios  (supersatura¬ 
tion  values),  the  barrier  becomes  finite  and  the  size  of  the  critical 
cluster  decreases.  At  a  loaarithmic  ratio  of  1.2,  a  barrier  no 
lonqer  exists,  and  at  higher  supersaturation  values  nucleation 
becomes  kinetically  controlled. 

This  is  consistent  with  a  purely  kinetic  formulation 
of  nucleation  which  was  developed  during  the  course  of  the  work. 

(See  Appendix  B,  Figure  2)  The  kinetics  of  cluster  breakup 
(evaporation,  kr)  are  plotted  along  with  the  forward  (clustering, 
kf)  rate  in  terms  of  its  product  with  the  concentration  of  the  con¬ 
densing  molecule  expressed  by  [^A-]  .  The  intersection  of  the  for¬ 
ward  and  reverse  curves  at  small  cluster  sizes  crives  rise  to  a 
stable  small  cluster  distribution.  At  [A-]*!,  corresponding  to 
100»  saturation,  no  curve  crossing  occurs  at  finite  cluster  sizes; 
nucleation  (curve  crossinq)  results  only  at  higher  (super)  satura¬ 
tions.  The  kinetically  controlled  regime  indicated  by  the  calcula¬ 
tions  shown  in  Figure  1  would  then  be  represented  by  the  case 
where  the  forward  curve  lies  above  the  reverse  one  in  Figure  2. 
Clearly,  a  complete  understanding  of  nucleation  requires  a 
consideration  of  both  reaction  kinetics  and  the  energy  barrier 
to  nucleation. 

Pertinent  to  the  problem  of  the  validity  of  the  Thomson 
equation  is  the  deqree  to  which  this  expression  reproduces  the 
requisite  thermodynamic  properties  of  the  ion  clusters.  Two 
systems  have  been  particularly  revealing  in  this  regard.  One 
concerns  the  clustering  of  water  to  the  alkali  metal  ions  and  Sr+, 


while  the  other  system  is  the  isoelectronic  molecule  ammonia 
clustered  to  similar  ions.  Utilizing  the  appropriate  derivatives 
of  the  Thomson  equation  we  have  compared  the  predicted  with  the 
experimental  thermodynamic  properties  deduced  from  ion  clustering 
studies . 

In  all  cases  the  Thomson  equation,  and  requisite  deriva¬ 
tives  thereof,  have  been  shown  to  well  represent  the  enthalpy  for 
successive  clustering  of  water  and  ammonia  to  ions  of  both  open 
and  closed  electronic  configuration.  Likewise,  considering  the 
entropies  for  the  clustering  of  water  molecules,  there  is  also  a 
reasonable  correspondence  between  the  values  predicted  by  the 
Thomson  equation  and  those  derived  experimentally  for  both  posi¬ 
tive  and  negative  ions.  (See  Figure  4  in  Appendix  D) .  However, 
a  consideration  of  the  entropy  for  ammonia  clustering  about  the 
alkali  metal  ions  reveals  a  qreat  discrepancy  between  theory  and 
experiment . 

Referring  to  Figure  6  (Appendix  B) ,  it  can  be  seen  that 
at  cluster  sizes  in  the  neighborhood  of  4  or  5  there  is  a  sig¬ 
nificant  departure  between  the  predicted  and  experimentally 
derived  entropy  values.  Not  only  is  the  departure  relatively  large 
at  this  deqree  of  clustering,  but  it  is  apparent  that  the  trend  is 
toward  an  even  greater  discrepancy  with  increasing  cluster  size. 
Clearly,  the  classical  formulations  do  not  hold  for  all  nucleating 
systems.  This  finding  is  important  in  that  it  reveals,  from  a 
molecular  point  of  view,  reasons  why  certain  systems  have  been 
found  to  be  reasonably  well  represented  by  the  energy  barrier 
predicted  by  the  Thomson  equation  while  others  display  consider¬ 
able  departure.  From  the  molecular  clustering  experiments  it  is 


now  apparent  that  this  discrepancy  is  in  fact  a  result  of  the 
nature  of  the  structure  of  the  prenucleation  embryos  to  ion  in¬ 
duced  nucleatior..  The  Thomson  equation,  by  its  very  nature,  im¬ 
plies  a  disordered  liquid-like  structure.  In  contrast,  the 
experimental  data  for  entropy  values  at  small  cluster  sizes  are 
seen  to  be  quite  negative,  becoming  more  so  with  successive 
clustering.  This  implies  that  the  small  ion  clusters  have  a 
relatively  ordered  structure,  which  provides  a  qualitative  picture 
of  the  reasons  why  the  Thomson  equation  is  not  valid  for  certain 
systems.  It  is  interesting  to  note  that  no  variation  in  the 
dielectric  constant  can  give  accord  between  the  predicted  and 
the  experimental  entropy  values.  But,  an  upward  adjustment  of 
the  surface  tension  for  the  ion  containing  droplets  with  decreas¬ 
ing  droplet  radius  was  found  to  give  somewhat  better  agreement. 

A  proper  formulation  of  surface  tension  for  a  droplet  containing 
an  ion  immersed  in  a  dielectric,  awaits  theoretical  attention. 

Consideration  has  also  been  given  to  the  utility  of  the 
data  in  interpreting  and  quantifying  the  interactions  between 
ion;  and  ligands  in  the  aqueous  phase  such  as  that  associated  with 
a  micro-droplet  aerosol  phase.  Details  are  given  in  Appendices 
D  and  G.  By  employing  the  high  pressure  mass  spectrometric 
technique  developed  in  our  laboratory,  we  have  been  able  to  study 
the  individual  interactions  at  selected  degrees  of  hydration. 
Thereby,  detailed  determination  of  bonding  and  the  mechanisms  and 
rates  of  reactions  can  be  obtained  for  ready  comparison  with 
theory.  Such  data  have  not  generally  been  attainable  from 
studies  typically  performed  in  bulk  solution,  and  their  utility 
in  interpreting  processes  involving  aerosols  has  been  limited.  As 
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detailed  in  Appendix  G,  a  particularly  important  result  of  the 
present  work  concerns  the  discovery  of  an  unexpectedly  simple  cor¬ 
relation  which  relates  cluster  bond  energy  to  bulk  phase  hydration 
energy,  as  a  function  of  cluster  size.  Despite  its  simple  form, 
it  has  been  found  to  be  valid  for  every  negative  ion  hydration 
reaction  for  which  data  are  available. 
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